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ABSTRACT 

An  error  probability  analysis  is  done  for  a  DFT  based,  M-ary  frequency-shift 
keying  (MFSK)  communications  system  employing  fast  frequency-hopped  spread 
spectrum  signals.  A  linear  combination  procedure  referred  to  as  noise -normalization  is 
employed  at  the  receiver  to  minimize  the  effects  of  partial-band  interference,  which  is 
modeled  as  additive  Gaussian  noise.  The  performance  of  the  receiver  is  studied  as  a 
function    of  signal   Doppler  shift  and  type  of  windowing    function    used  in  the  DFT. 

The  use  of  fast  frequency-hopped  spread  spectrum  is  found  to  improve  the 
performance  of  the  DFT  based  receiver  in  all  but  the  most  severe  cases  of  Doppler  shift. 
The  use  of  a  non-uniform  window  (i.e.,  a  Hamming  window)  to  improve  receiver 
performance  is  effective  only  in  the  presence  of  large  Doppler  shifts.  The  amount  of 
Doppler  shift  necessary  to  warrant  the  use  of  a  non-uniform  window  depends  on  the 
amount  of  jamming  noise  power  at  the  receiver,  but  is  relatively  insensitive  to  the 
frequency-hop  rate  used.  In  general,  in  the  absence  of  any  information  concerning  the 
nature  of  the  received  signal,  a  non-uniform  window  should  be  used  because  the 
performance  degradation  experienced  at  small  Doppler  shifts  is  insignificant  compared 
to  the  performance    enhancement    gained   at  large  Doppler  shifts. 
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I .   INTRODUCTION 

A.   NON-UNIFORM  WINDOWING  OF  THE  DISCRETE  FOURIER  TRANSFORM 

The  use  of  the  discrete  Fourier  transform  (DFT)  to 
noncoherently  detect  signals  employing  M-ary  frequency-shift 
keying  (MFSK)  modulation  has  gained  popularity  recently  due  to 
advances  made  in  high  speed  digital  systems.  At  the  current 
level  of  technology,  however,  the  utility  of  the  DFT  based 
receiver  remains  limited  due  to  the  speed  of  real-time  DFT 
algorithms.  Even  the  performance  enhancement  obtained  by  using 
the  fast  Fourier  transform  (FFT)  is  not  sufficient  to  overcome 
this  limitation.  However,  certain  applications  -  namely  fast 
frequency-hopped  spread  spectrum  systems  -  can  effectively 
utilize  the  DFT  receiver  because  lower  data  rates  frequently 
correspond  to  higher  immunity  to  narrowband  interference. 
[Ref.  l:p.  2014] 

Problems  arise  in  using  the  DFT  based  receiver  if  the 
transmitted  frequency  does  not  match  the  assigned  frequency  of 
the  DFT  bin.  This  situation  results  when  the  signal  frequency 
is  different  than  that  required  by  the  detector,  as  when  the 
signal  experiences  a  Doppler  shift  prior  to  detection.  This 
leads  to  a  situation  where  transmitted  signal  energy  falsely 
contributes  to  an  incorrect  DFT  frequency  sample:  a  phenomena 
known  as  self-induced  crosstalk  [Ref.  l:p.  2015].  This  problem 


is  not  specific  to  the  DFT  based  receiver.  What  is  unique 
about  this  type  of  receiver  is  its  ability  to  minimize  self- 
induced  crosstalk  through  the  use  of  non-uniform  windows.  An 
example  of  a  DFT  based  detector  experiencing  self-induced 
crosstalk  is  presented  in  Fig.  1.1.  Figure  1.1(a)  depicts  the 
situation  where  the  received  frequency  is  perfectly  matched  to 
the  first  output  bin  of  the  detector ,  while  Fig.  1.1(b)  shows 
the  results  of  slightly  modifying  the  received  signal's 
frequency.  As  can  be  seen,  the  magnitude  of  the  bin 
corresponding  to  the  signal  is  reduced,  while  the  magnitudes 
of  neighboring  bins  which  represent  other  frequencies  are 
significantly  larger.  Large  amounts  of  frequency  deviation  can 
easily  lead  to  unacceptable  bit  error  rates. 

One  method  that  is  used  to  minimize  the  effects  of  self- 
induced  crosstalk  involves  the  use  of  nonuniform  windowing; 
that  is,  time-sampled  data  is  not  equally  weighted.  The 
frequency  response  of  typical  nonuniform  windows  have  broader 
main  lobes  than  do  rectangular  windows  (uniform  weighting)  but 
exhibit  sharper  rolloff  characteristics.  The  sharp  rolloff 
characteristics  of  nonuniform  windows  partially  compensate  for 
the  effects  of  self-induced  crosstalk,  but  in  some  cases  the 
existence  of  a  broader  main  lobe  creates  distortion  where  it 
would  not  otherwise  exist  [Ref.  l:p.  2015].  Nonetheless, 
nonuniform  windowing  can  be  used  successfully  to  combat  large 
doppler  shifts  in  low  data  rate  environments  as  demonstrated 
in  [1]  and  [2]. 
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Figure  1.1:  Outputs  of  a  DFT  (a)  without  a  doppler  shifted 
signal,  and  (b)  with  a  doppler  shifted  signal.  The  adverse 
effects  of  a  doppler  shifted  signal  can  severely  degrade 
communications . 


B.   FACTORS  AFFECTING  RECEIVER  PERFORMANCE 

1.   Fast  Frequency-Hopping  and  Partial-Band  Interference 

When  the  carrier  of  the  transmitted  frequency  is 
changed  in  an  apparently  random  manner  during  the  transmission 
of  a  signal,  the  overall  effect  is  to  spread  the  energy  of  the 
signal  over  a  large  bandwidth.  This  spreading  technique  is 
called  frequency-hopping  and  is  illustrated  in  Fig.  1.2.  Two 
types  of  frequency-hopping  are  considered:  fast  and  slow 
frequency-hopping.  Slow  frequency-hopping  systems  change,  or 


hop,  the  carrier  once  per  data  symbol,  while  fast  frequency- 
hopping  systems  hop  the  carrier  more  than  once  per  data 
symbol.  The  higher  the  hopping  rate,  the  more  the  frequency 
spectrum  approximates  that  of  true  spread  spectrum  systems 
utilizing  pseudonoise  (PN)  codes  [Ref.  3:  pp.  10-14]. 
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Figure  1.2:  An  example  of  hopping  the  carrier  frequency  of 
a  signal. 


The  intent  of  a  hostile  jammer  is  to  disrupt 
communications.  In  order  to  effectively  jam  a  frequency-hopped 
signal,  a  smart  jammer  may  employ  partial-band  interference; 
although,  this  type  of  interference  can  also  be  caused  by 
other  unintended  narrowband  sources.  In  utilizing  partial-band 
interference,  the  jammer  concentrates  all  of  the  jamming 


energy  into  a  fraction  of  the  total  bandwidth  and  randomly 
hops  this  signal  through  the  spectrum  in  the  same  fashion  as 
the  signal  employing  frequency-hopping.  By  doing  this,  the 
frequency  hopped  signal  cannot  always  avoid  the  portion  of  the 
spectrum  that  is  jammed;  the  jammer  must  be  dealt  with  via  the 
receiver. 

The  assumptions  involved  in  the  analysis  to  follow  are 
that  the  jammer  is  limited  in  power;  that  the  interference  can 
be  modeled  as  additive  Gaussian  noise;  and  that  the 
interference,  when  present,  is  present  in  all  detection 
channels  with  probability  y .  Hence,  the  probability  that 
interference  is  not  present  in  any  branch  is  1-y.  Thus,  y 
represents  the  fraction  of  the  spread  bandwidth  that  is 
jammed.  If  the  partial-band  interference  is  assumed  to  have  a 
power  spectral  density  of  N./2  over  the  entire  spread 
bandwidth,  the  amount  of  interference  present  at  the  receiver 
due  to  partial-band  interference  is  y~1N./2  with  probability  y 
and  0  with  probability  1-y.    [Ref.  4:p.  3] 

In  addition  to  partial-band  interference,  the  system  is 
assumed  to  be  corrupted  by  additive  white  Gaussian  noise  with 
a  power  spectral  density  of  NQ/2 .  Accordingly,  the  power 
spectral  density  of  the  total  noise  at  the  receiver  input  is 
N0/2  +  y"1N./2  with  probability  y  and  NQ/2  with  probability  1-y. 
If  the  receiver  equivalent  noise  bandwidth  is  W  Hz,  then  the 
total  noise  power  at  the  receiver  is  NQW  with  probability  1-y 
and  is  (  y~1N.+N0)W  with  probability  y.     [Ref.  4:p.  3] 


2 .   Fading 

In  addition  to  partial-band  interference,  the 
signal  is  assumed  to  be  further  degraded  by  fading.  Fading  is 
caused  by  components  of  the  signal  arriving  over  different 
paths  which  combine  to  form  the  total  detected  signal.  Due  to 
the  different  path  lengths  traversed  by  the  various  signal 
components,  the  combined  signal  amplitude  may  be  greater  than 
or  less  than  the  expected  amplitude  without  fading  due  to 
constructive  and  destructive  interference.  Fading  severely 
degrades  the  performance  of  a  real  world  communications  system 
and  must  be  taken  into  account  in  any  plausible  receiver 
design  and  analysis. 

For  this  analysis,  the  signal  is  assumed  to  experience 
slow,  frequency  non-selective,  Rician  fading.  Slow  fading 
refers  to  the  fact  that  the  amplitude  of  the  signal  is  assumed 
to  remain  constant  at  least  over  the  duration  of  one  hop, 
while  frequency  non-selective  implies  that  all  frequency 
components  of  the  signal  experience  identical  fading.  Rician 
fading  is  characterized  as  having  a  diffuse  signal  component 
and  a  direct  signal  component.  Assumme  that  the  amplitude  of 
the  received  signal  in  a  Rician  fading  channel  is  a/2.  Then 
the  probabilty  density  function  of  the  random  variable  A  is 
given  by 


A_e  >°f2    T(J*L)u(a).  <D 


fA(a)=^-e      <"*     l0(-^-)u(a)t 


where  u(-)  is  the  unit  step  function,  r2  is  the  power  in  the 
direct  component  of  the  signal,  2af2  is  the  power  in  the 
diffuse  component,  and  I0(  •  )  is  the  modified  Bessel  function 
of  order  zero.  As  can  be  seen,  (1)  reduces  to  the  well  known 
Rayleigh  distribution  when  the  value  of  the  direct  component 
is  zero.  [Ref.  5] 

C.   THE  NOISE -NORMALIZED  DETECTOR 

A  noise-normalized  detector  is  a  linear  receiver  that  uses 
the  received  noise  power  to  normalize  the  output  of  each 
branch  before  a  decision  is  made  regarding  which  signal  is 
present.  This  type  of  receiver  is  combined  with  fast  frequency 
hopping  to  improve  overall  system  immunity  to  partial-band 
interference.  The  noise-normalization  is  implemented  before 
the  hop  receptions  are  combined  to  form  the  detection 
statistics.  This  type  of  detection  scheme  tends  to  lessen  the 
influence  on  the  overall  decision  statistic  of  each  hop  when 
interference  is  present,  while  increasing  the  influence  of 
each  hop  when  no  interference  exists  [Ref.  4:p.  2].  Because 
the  output  statistics  of  each  branch  of  the  receiver  are 
independent  of  one  another,  closed  form  solutions  for  the 
receiver  performance  can  be  obtained. 

This  thesis  studies  the  effects  of  nonuniform 
windowing  on  a  noise-normalized  DFT  based  receiver  where  the 
signal  experiences  partial-band  interference,  Doppler  shift, 
and  Rician  fading.  The  effects  of  fast-frequency  hopping  the 


carrier  frequency  are  analyzed  and  the  results  for  various  hop 
rates  discussed.  As  in  [1]  and  [2],  a  statistical  analysis  is 
done  to  determine  the  amount  of  frequency  deviation  necessary 
to  warrant  the  use  of  a  nonuniform  window  over  a  rectangular 
window. 


II.  ANALYSIS  OF  THE  NOISE-NORMALIZED  DFT  BASED  DETECTOR 

The  receiver  under  consideration  is  shown  in  Fig.  2.1.  The 
incoming  signal  hop  is  first  brought  back  to  baseband  and  then 
sampled  at  the  Nyquist  rate.  The  DFT  of  the  signal  is  computed 
using  the  appropriate  windowing  function  and  the  magnitude  of 
each  output  bin  is  squared.  Hence,  in  the  absence  of  both  a 
Doppler  shift  and  a  nonuniform  window,  this  receiver  is 
functionally  equivalent  to  a  conventional  MFSK  receiver 
employing  quadratic  detectors.  At  the  same  time,  the  total 
noise  power  at  the  receiver  is  measured,  modified 
appropriately  by  the  windowing  function  of  choice,  and 
inverted  to  form  the  noise-normalization  term.  Each  bin  output 
is  multiplied  by  the  noise-normalization  term.  The  decision 
statistics  are  then  obtained  by  summing  over  all  hops 
comprising  the  transmitted  symbol.  Finally,  the  detection 
statistics  are  compared  and  the  largest  one  used  to  determine 
the  estimate  of  the  symbol  sent.  Errors  result  when  the 
detection  statistic  in  a  channel  with  no  signal  is  greater 
than  the  detection  statistic  of  the  signal  channel.  A  detailed 
analysis  of  each  part  of  this  design  follows. 


Figure  2.1:  The  DFT  based  noise-normalized  receiver 
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A.   DETERMINATION  OF  THE  PROBABILITY  DENSITY  FUNCTION  FOR  EACH 
DFT  OUTPUT  BIN 

The  signal  input  to  the  receiver  in  Fig.  2.1  every  Th 
seconds  is  of  the  form 

r  ( t)  =av/2cos  [2n  (fe+fm+f)  t]  +n'{t)  ,  (2) 

where  CKt£Th,  a  is  a  Rician  random  variable  having  the 
probability  density  function  of  (1),  a/7  is  the  Rician 
distributed  signal  amplitude,  fc  is  the  carrier  freguency 
hopped  every  Th  seconds,  fm  is  the  intended  symbol  freguency, 
f  is  a  Doppler  shift,  and  n'(t)  is  zero  mean,  additive 
Gaussian  noise.  The  signal  freguencies  are  chosen  to  be 
orthogonal  during  the  time  interval  0<_t<Th  and  are  obtained 
using  the  relation 

fa>=—  + .  (3) 

1h  1h 

where  m'=l,2,..,M  and  Af  is  an  integer  representing  the 
spacing  between  each  of  the  M  signals.  To  simplify  the  signal 
representation,  the  following  substitutions  are  made: 

gAl+Af(/n/-l)  (4) 

and 

ekf'Th  (5) 

so  that  after  the  low-pass  filter  (2)  becomes 
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r'it)  =a/2cos[  2n  {q+e)  l]  +n(t)  ,  (6) 


where  e  is  the  fraction  of  Doppler  shift  with  respect  to  the 
bandwidth  of  the  baseband  information  signal,  q  is  an  integer 
related  to  the  signal  frequency  f  , ,  and  n(t)  is  the  filtered 
noise.  The  discrete  form  of  the  received  signal  is  obtained  by 
using  a  sampling  rate  of  N/Th,  where  N  is  the  number  of 
samples  taken  during  the  symbol  reception: 

r(p)  =av/2cos  [  2*<g+e>P]  +n(p)  ,  (7) 

N 

where  p=0,l,..,N-l  and  is  the  discrete  signal  sample.  [Ref. 
2:pp.  7-9] 

The  output  of  each  bin  is  a  complex  quantity  and  can  be 
separated  into  real  and  imaginary  parts: 

Yjk\a)  =r/aB(Jc:|a)  +jYIm(k\a)  ,  (8) 

where  k,  which  ranges  from  0  to  N-l ,  is  an  integer 
representing  the  DFT  output  bin  and  Yn(k\a) ,  ^(Ala),  and 
Ylm(k\a)  are  respectively  the  DFT  output,  the  real  part  of  the 
DFT  output,  and  the  imaginary  part  of  the  DFT  output  for 
channel  k  and  frequency  hop  m  conditioned  on  the  Rician 
distributed  random  variable  a.  The  real  and  imaginary  parts  of 
the  DFT  are  calculated  separately 
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YSa(k\a)=Y^10r(p)w(p)cos[^P]  (9) 

=£^  [v^acos  [  27r  {q:*)p]  +n(p)  ]  w(p)  cos  [^£]  , 

and 

li,(*|a)-I^r(x.)v(p)Bln[^B]  (10) 

=E£  [v/2acos-[  2"  'g+e»P]  +n(p)  ]  „(P)  sin[.g^E]  , 

where  w(p)  represents  the  windowing  function.  Assuming  slow 
Rician  fading  so  that  the  amplitude  of  the  signal  remains 
constant  during  each  hop  interval,  the  conditional  probability 
density  function  for  the  random  variable  Xkm  is 

-Km  \  (11) 

rxkm\A^Xkm\a) "€ 

where  Xkm  is  the  magnitude  of  bin  k   and  hop  m   given  by 

X,   =\Y  (k\a)  I 

(12) 


-(-xw'-nL) 

o  2 

X 

2ox2          T    (   Xkm\Xkm 
O 

X 

X*n=\Yjk\a) 


jukm  is  the  magnitude  of  the  expected  value  of  bin  k     and 
frequency  hop  m   given  by 


and  a  2  is  the  variance  of  bin  k   given  by 


(13) 


ax2  =  ±E[\Yn(k\a)  -Yjk\a)\2]  .  (14) 


[Ref.    l:p.    2016] 
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The  mean  of  Xkm  for  a  fixed  amplitude  a  is  easily- 
determined  by  observing  that  the  only  random  quantity  is  n(p) 
which  is  zero-mean.  This  leads  to  the  conclusion  that  the  mean 
is  simply  the  output  of  the  DFT  without  noise.  The  real  and 
imaginary  parts  of  the  mean  conditioned  on  a   are 


yxm~*#EZ*ip)coBi  2*^+e)  >  cos  <  »5e*) 


(15) 


™Rm 


and 


^^w-^jt:  »<*>  c°s  <  2itpire)  > sin  <  -^ > 


'P=o  tf  itf    .    (16) 


*%m 


An  expression  for  /Jkm2  conditioned  on  a  is  then  determined  by 
substituting  (15)  and  (16)  into  (13)  and  squaring  both  sides: 

=a2(m'2R+m>2I)  (17) 

where  /3km  is  the  magnitude  squared  of  the  mean  of  the  DFT 
output  without  fading  when  the  signal  power  is  normalized  to 
unity,  that  is, 

P-=2  <  E£  "<P>  cos  <  27tP^+6)  )  cos  ( J5E* )  ]  a 

+  [Vv-1cos(27t^((3r+e)sin(^P^)]2). 

[Ref.  2:pp.  10-11] 

The  variance  of  Xkm  is  obtained  by  substituting  (9),  (10), 
(15),  and  (16)  into  (14)  and  simplifying 
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ax2=^E[|Fi2m(ic|a)  -Y^ikla)  +j{YIm{k\a)  -Ylw(k\a)  )  |2] 

=E[\yy-l  w(p)  n  (p)  cos  (  i^2  )  ( 19 ) 

+-7'£^  "<p>  *  (p)  sin  < -^^  H2]  • 


The  magnitude  squared  of  (19)  becomes 


+(C,(p),(P)Sin(i|P))2]  (20) 

~BlY£tt"W  "I*)  eoB(^)coei^)  nip)  nkB) 
+T£J^*lP)"ls)Bin(^)eiiL{*2&)n{p)n(s)], 


(21) 


where  s  is  an  integer  having  the  same  range  of  values  as  p. 
Because  the  expected  value  operator  is  a  linear  operator  and 
the  only  random  quantities  are  the  noise  functions,  (20) 
simplifies  to 

o^ly^'lj^-1  w(p)  w(s)  cos  (^-) 

x       2  ^— 'p=o  a-^s=o  pj 

•  cos{H£Z)  E[n{p)n(s)] 

N 

+±Y"-lY"-y(p)  wis)  sin(2*£P) 

2  ^— 'p=o  ^^s=o  j\f 

•  sin(-^MjF[n(p)73(fir)]  , 

N 

The  autocorrelation  function  is  defined  as 

R^ip^)  =E[n(p)n(s)]  ,  (22) 

which,  because  the  noise  functions  represent  white  noise,  is 

R^ip.s)  =WNtb(p-s)  ,  (23) 

where  6(  ■  )  is  the  Dirac  delta  function,  W  represents  the 
equivalent  noise  bandwidth  of  the  low-pass  filter  at  the  front 
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of  the  receiver,  and  Nt  is  the  power  spectral  density  of  the 
total  noise  seen  at  the  receiver  [Ref.  4].  The  quantity  Nt 
equals  either  NQ  when  jamming  is  not  present  or  NQ+N./y  when 
jamming  is  present.  The  expected  value  needed  to  solve  (21)  is 
obtained  by  substituting  (23)  into  (22): 

E[n(p)n(s)]=WNt6(p-s)  .  (24) 

The  final  expression  for  the  variance ' of  Xkm  is  obtained 
by  substituting  (24)  into  (21)  and  simplifying: 

«,'--^I^ES-Cp)w(«)cob(^E)cos(^)»(p-) 

l^trv-i  v^-i  „/~n  „/„\  „.:~  /  2nkps  „.     ,  2nks 


+J^ty^iy^iw{  p)w{s)3±n{2jzkp)s±n{2nks)b{       , 

2  L^ip=Q  Z-js=q       ^  N  N  (25) 

£*-;  ^  (P)  [cos^  ( ^£p )  +sin2  ( .^p )  ] 

*— 'E>=0  A/  A; 


WATt 


2  ^p=°  iV  AT 

2™£r"*«(p>. 


'p=0 

The  variance  in  each  bin  is  affected  by  both  the  power 
spectral  density  of  the  external  noise  and  the  choice  of 
windowing  function  used  in  the  DFT. 

The  probability  density  function  for  the  output  of  each 
bin  is  found  by  substituting  (17)  into  (11)  and  integrating 
over  the  range  of  the  Rician  distributed  amplitude: 
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_  C"  xkm  2ox2  j    /   axknflJPkm\ 


x 

-(a^r2) 


e  Jo< — 2)da  (26) 


of2  of 

Xkm  2      0jt2      0f2 


O    2C f2 
x      r 


km    .        1 


_32  <      fkm    +       J.       .  

■  fae        ><    -/  x0(axtai^)X0(^)da 
Jo  ao* 

x  f 


By  making  the  substitutions  I0(x)=J0(  jx) ,  and  by  defining 


$km     ,         1 


0^4        —    +■ 


2o  z      2o/ 

«     "      2  2  <27) 

Pj™°f    +°x 


2o/of2 


_<■!«£, 


and  using  the  integral  (from  [7]) 

rxe-^\jv(ax)Jv(£x)c^=^e~  ~^'lv(^),     (28> 
Jo  202  2£)2 

the  probability  density  function  is 


f  „  (x,J  =  ■ 


■X^^/on 


XJcm  2ox2or2 


2^<V°/  (29) 

1      .   (    xkm$kir,  ^     T2    }  


2Qdox'of 

The  next  step  involves  substituting  (27)  into  (29)  which  gives 
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*j>/2-r2°x2  ^     «&P*d»/+r»o,« 


f        (x      )  -  (  X*m  )  e         2ox2o/  2ox2o/(pitmo/+ox2) 


*)cra 


Equation  (30)  simplifies  to 


$knPf2+°x  (30) 


__1   ,     xton4'PfaB^    j 

-XVm  2 


^a e       i^M'  j  (   * taT/P^  )  ,         (31) 


In  order  to  clarify  the  notation,  the  following  substitutions 
are  made: 

a^^of+o*.  (32) 

and 

«^r/^.  (33) 

With   these   substitutions,   (31)   is   seen   to   be   Rician 
distributed 


-^km2*akm2) 

-*Jan         Km           n       2 
°Jcm 

2o  te2 

T    (  Xkmakm 

0                    2 
°km 

), 


(34) 


and  reduces  to  the  result  obtained  in  [1]  for  the  case  of  a 
nonfading  channel  (2af2=0  and  r=l).  [Ref.  2:pp.  11-12] 
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B.   ANALYSIS  OF  THE  NOISE-NORMALIZED  RECEIVER 

1.   Determination  of  the  Probability  Density  Function  for 
the  Detection  Statistic 

As  illustrated  in  Fig.  2.1,  the  random  variable  Zkm   is 
given  by 

Z*n=^4.  (35) 

0 

X 

where  Zkm  is  the  random  variable  for  bin  k  and  frequency  hop 
m.  Because  all  the  random  variables  Xkm  are  independent,  the 
transformation  of  random  variables  is  accomplished  by  using 
the  relation  (from  [8]) 


*3   (Xkm(0)  ' 


**.<»«■>-    r.     '.  -  <36> 


where 

Ztm'ffiXtm).  <37> 

g'  (xkm)  is  the  derivative  of  g(xkm)  with  respext  to  xkm,  and 
xkm(0)  ^s  t^e  real  root  of  (35).  The  denominator  in  (36)  is 
found  by  substituting  (35)  into  (37)  and  differentiating: 
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a'(x         \=.2x**><°> 


O   2 
x 


=2 


Zkm 


(38) 


\  o 


2 

X 


The  transformation  of  random  variables  for  the  noise- 
normalized  detection  statistic  is  determined  by  substituting 
(38)  into  (36): 


fZka^Zkm'       2  \|  z 


£^<i/57*j  •         (39) 


The  probability  density  function  of  Zkm     is  obtained  by 
substituting  (34)  into  (39): 


(°x2zj™+«L> 


f      [y     )-    Q*     r  2o^         T   (  a*»° */***)  (40) 

The  actual  detection  statistic  needed  for  the 
comparison  is  Zk,    which  is  defined  as 


vZL^  <41> 


where  L  is  the  number  of  freguency  hops  per  symbol.  Since  some 
hops  are  jammed  while  others  are  not,  the  probability  density 
function  of  Zk  conditioned  on  the  number  of  jammed  hops  is 
simply  the  multiple  convolution  of  the  probability  density 
functions  of  the  jammed  and  unjammed  hops.  For  instance,  if  a 
symbol  has  three  hops  and  two  of  them  are  jammed,  the 
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probability  density  function  of  the  detection  statistic  is 
determined  as  follows: 

fZk(zk)=fzJzkl)J™^fzJzk2)J™e*®fzJzk3)^*™*d,        (42) 

where  0  represents  convolution.  In  general,  the  probability 
density  function  of  Zk  conditioned  on  I  out  of  L  hops  jammed 
is 

fZk  (**>-[  fZkm  < **■>  ^^  ®  J®  [ ^  ( z**)  unjammed]  ®  <L~J) .        < 43 ) 

where  ®cn  represents  a  cn-fold  convolution.  [Ref.  4:p.  7] 
2.   Determining  the  Probability  of  Bit  Error 

Once  the  random  variable  transformation  is 
accomplished,  the  next  step  involves  calculating  the 
probability  of  bit  error.  An  error  will  occur  if  the  power  in 
a  channel  other  than  the  signal  channel  is  greater  than  the 
power  in  the  signal  channel.  Mathematically,  this  is  expressed 
as 

PriZ^Zj)  =fXfz.(zj)fZifZi(zi)dzidzjl  (44) 

where  i,j=l,...,M  and  i*j;  and  the  signal  is  assumed  to  be  in 
branch  i.[Ref.  4:  p.  7].  The  double  integral  in  (44) 
determines  the  probability  that  the  power  present  in  the 
signal  channel  is  less  than  the  power  present  in  a  noise 
channel . 
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Because  the  distribution  of  power  in  each  noise 
channel  is  not  identical  when  Doppler  shift  is  present,  the 
comparison  in  (44)  must  be  carried  out  for  each  noise  channel 
and  the  results  combined  to  give  the  total  probability  of 
error.  An  upper  bound  to  the  total  conditional  probability  of 
error  is  found  using  the  union  bound: 

PSn^J-iPriZ^Z^  .  (45) 

where  Psn  is  the  probability  of  symbol  error  given  1  hops  are 
jammed  with  the  signal  present  in  channel  i.  [Ref.  l:p.  2017] 
Because  the  distribution  of  power  in  each  bin  varies 
with  the  location  of  the  signal,  the  amount  of  Doppler,  and 
the  relative  direction  of  the  freguency  shift,  an  average  must 
be  taken  over  all  signal  channels  and  direction  of  Doppler 
shift: 

Ps2  (e)  =-^5Xi  Psn  (  +«)  +Psn  (-e)  ,  (46) 

where  Psl(e)  is  the  total  probability  of  symbol  error  given  1 
hops  are  jammed  conditioned  on  the  signal  experiencing  a 
Doppler  shift  of  e  and  Psn(-)  is  the  probability  of  symbol 
error  given  1  hops  are  jammed  with  the  signal  present  in 
channel  i  conditioned  on  the  signal  experiencing  a  Doppler 
shift  of  (+£)  or  (-£).  The  union  bound  on  the  conditional 
probability  of  symbol  error  is  calculated  by  combining  (45) 
and  (46): 
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pSi(e)<^-£"  £*=1  [Pr(Zi<ZJ-|+e)+Pr(Zi<ZJ|-G)]  .    (47) 

The  total  conditional  probability  of  symbol  error  assuming  1 
hops  are  jammed  is  calculated  by  combining  (44)  with  (47): 

**!<«> '-sE-aEja  f//z,i+e(^i^)/o%l+e(zii+G)dzidz,  (48) 

+/o%l-  ( z*  I  "€)  /o  %N ( zi  I  "€)  dzidzi]  ■ 

[Ref.    l:p.    2017] 

The  total  probability  of  symbol  error  for  the  noise- 
normalizing  receiver  in  the  presence  of  partial-band 
interference  is 


Ps(e)=5^=o|^)YJ(l-Y)L-J^J(€)/  (49) 


where  Ps{e.)  is  the  total  probability  of  symbol  error  [Ref. 
4:p.  6].  To  determine  the  worst  case  probability  of  symbol 
error,  (49)  is  evaluated  for  different  values  of  y,  the 
results  compared,  and  the  highest  conditional  probability  of 
symbol  error  recorded.  The  worst  case  probability  of  bit  error 
is : 

Pb(e)=-^Ps(e)  ,  (50) 

M-l 

where  Pb(e)    is  the  probability  of  bit  error. 

An  exact  solution  for  the  conditional  probability  of 
symbol  error  is  calculated  by  evaluating  the  probability  of 
not  making  an  error  (i.e.,  the  probability  that  the  power  in 
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all  noise  channels  is  less  than  the  power  in  a  signal 
channel) : 

PziZ^Z,)  =f"fZi(zi)foZifZj(zj)dzjdzi.  (51) 

The  total  conditional  probability  of  not  making  an  error  is 
found  by  the  following  relation: 

Psn(e)=llUPr(Zj<Zi\e) 

J**  (52) 

=n>i  Cf  *  iz^e)  L'% (  zj  ie)  dzidzi ' 


where  Ps-l(z)  is  the  probability  of  not  making  an  error  given 
that  the  signal  in  channel  i  with  J  hops  jammed  and  the  signal 
experiencing  a  Doppler  shift  of  e.  The  probability  of  making 
an  error  given  that  the  signal  is  in  channel  i  with  I  hops 
jammed  is  then  given  by 

Psn(e)  =1-Psu(e)  .  (53) 

which  is  then  used  in  (46),  (47),  (48),  (49),  and  (50)  to 
obtain  the  exact  probability  of  bit  error.  Due  to  the  enormous 
amount  of  computing  power  needed  to  numerically  integrate  (51) 
and  the  desire  to  compare  the  results  with  those  obtained  in 
[1]  and  [2],  the  conditional  probability  of  error  based  on  the 
union  bound  was  chosen  for  this  thesis. 
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III.   NUMERICAL  PROCEDURE 

A.   PARAMETERIZATION  OF  VARIABLES 

In  order  to  compare  the  performance  of  the  noise- 
normalizing  receiver  with  a  conventional  DFT/MFSK  receiver, 
the  number  of  points  used  in  the  DFT,  the  cut-off  frequency  of 
the  low  pass  filter,  and  the  spacing  between  the  M  signals  are 
chosen  to  be  consistent  with  those  used  in  [1]  and  [2].  The 
following  definitions  result  from  adopting  these  conventions: 

Nk4M,  ( 54 ) 

W&2M,  (55) 

and 

Af±2,  (56) 

where  N  is  the  number  of  signal  samples  per  hop,  M  is  the 
modulation  order,  W  is  the  effective  noise  bandwidth,  and  Af 
is  an  integer  representing  the  message  frequency  spacing. 

A  hop  time,  Th,  of  one  second  results  from  sampling  the 
received  signal  at  the  Nyquist  rate.  When  this  value  of  Th  and 
(56)  are  substituted  into  (3),  it  is  determined  that  fm 
consists  only  of  odd  values  ranging  from  1  to  2M-1 . 
Furthermore,  the  spacing  of  frequencies  in  the  DFT  is 
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N 
_  4M  (57) 

4M 

where  4f  is  the  frequency  resolution  of  the  DFT  [Ref.  9]. 
Because  of  these  choices  one  empty  bin  exists  between  each 
possible  DFT  signal  output. 

A  parameter  of  concern  is  the  energy-per-bit  to  thermal 
noise  power  spectral  density  ratio  defined  as  Eh/NQ,  where  Eb 
is  the  transmitted  energy  per  bit.  The  signal-to-noise  ratio, 
which  is  a  measure  of  the  ratio  of  signal  power  to  noise 
power,  is  defined  as 


SNRk 


r2+2o  2 

1  Z°f    .  (58) 


N0W 


For  an  M-ary  system,  symbols  are  sent  instead  of  bits.  The 
power  per  symbol  is  the  sum  of  the  power  in  the  direct 
component  and  the  power  in  the  diffuse  component: 

Es=(2of2+T2)Ts,  (59) 

where  Es  is  the  transmitted  energy  per  symbol  and  Ts  is  the 
symbol  duration  time  related  to  the  hop  duration  time  by 

Ts=LTh.  (60) 

The  energy-per- symbol  is  related  to  the  energy-per-bit  by 
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Es=Eblog2  (M)   .  (61) 

All  the  variables  used  in  this  thesis  are  expressed  as  per-bit 
quantities;  therefore,  a  relationship  between  the  energy-per- 
hop  and  the  energy-per-bit  is  needed.  Because  each  symbol  is 
comprised  of  L  bits, 

Eh=-5s,  (62) 

h      L 

where  Eh  is  the  energy-per-hop.  A  relationship  between  the 
energy-per-bit  to  noise  power  spectral  density  ratio  and  the 
signal-to-noise  ratio  is  obtained  by  combining  (58)  with  (59) 

— £=- ?—SNR.  (63) 

N0      log2  (M) 

Proper  choice  of  the  symbol  duration  time  (i.e.,  Ts=l)  and  the 
equivalent  noise  bandwidth  (i.e.,  W=log2(M)  )  gives  Eb/NQ=SNR. 
In  general,  however,  these  two  ratios  are  not  equal  and  must 
be  kept  separate  in  the  receiver  analysis. 

An  expression  relating  the  energy-per-bit  to  noise  power 
spectral  density  ratio  to  the  signal-to-noise  ratio  is 
obtained  by  combining  (58),  (60),  and  (63)  and  simplifying: 

Eb_   (W)   (L)  (Th)    r2+2of2  (64) 

N0  log2(M)  N0W 

For  the  case  of  a  receiver  operating  in  a  jammed  environment, 
NQ  is  replaced  with  NQ+N-/y  when  the  hop  is  jammed  and  (64) 
becomes 
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Eb{  1      ,    r2+2c/    (^)(L)(TA) 

iV0  l+Nj/ (yN0)         WiNo+Nj/y)       log2(M) 

The  relationship  between  the  energy-per-bit  to  noise  power 
spectral  density  and  the  signal-to-noise  ratio  is  a  function 
of  the  effective  noise  bandwidth  of  the  system,  the  number  of 
freguency  hops  per  symbol,  the  hop  duration  time,  and  the 
level  of  M-ary  signalling  used. 

Another  parameter  for  analysis  is  the  direct  to  diffuse 
ratio 


X2 

2a 


DTDk-^.  (66) 


f 

The  following  expressions  for  the  direct  and  the  diffuse 
signal  power  components  result  from  combining  (65)  and  (66) 
and  simplifying: 

af2  = ^ (67) 

f      2(l+DTD)Th 


and 


T2  =  E±{     DTD     }  (68) 

Th      1+DTD 


In  doing  this  analysis,  the  values  of  Eb/NQt  N-/NQ/  y,  W,  Th, 
and  DTD  are  assumed  and  used  to  determine  the  values  of  r2  and 
cf*. 


28 


B.   CONVOLUTING  THE  PROBABILITY  DENSITY  FUNCTIONS 

In  the  case  of  fast  frequency-hopping,  the  probability 
density  function  of  the  total  detection  statistic  is  the  L- 
fold  convolution  of  the  probability  density  functions  of  the 
detection  statistic  of  each  hop.  As  shown  in  (49),  the  total 
probability  of  error  represents  all  possible  combinations  of 
jammed  and  un jammed  hops.  Based  on  this,  a  total  of  2M(L+1)  L- 
fold  convolutions  are  required  for  each  possible  signal  branch 
and  for  each  combination  of  jamming  power  and  fraction  of 
bandwidth  jammed  that  are  studied.  These  convolutions  are  not 
required  in  the  case  of  slow  hopping  because  the  total 
detection  statistic  is  equal  to  the  detection  statistic  of  the 
single  hop.  These  multiple  convolutions  are  generated  by 
taking  advantage  of  the  properties  of  the  Laplace  transform. 

The  Laplace  transform  of  the  probability  density  function 
of  Z,_   is 


lm 


g[fzln)(Z^))]=Fz(n)(s) 
*-im  ''im 

-fmfxj*lz}?)e—Pdz} 

Jo         lm 


<69> 


where  n=l,2  and  n  =  l  represents  the  condition  of  a  jammed  hop 
and  n=2  represents  the  condition  of  an  unjammed  hop.  The 
Laplace  transform  of  interest  is  obtained  by  substituting  (40) 
into  (69) 
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(o  2z(n)+a.  2<n)) 


Jo  °^     „*to)  }       iw 


Equation  (70)  is  evaluated  by  defining 

7<x  •  (n)  a 
B=   "    "        .  (71) 

and 


o  2 


C2  =  s+  — f—  ,  (72) 


->„2(j}) 


and  making  the  substitutions  u2=z-m   and  I0(x)=JQ(  jx) ,  and  using 
the  integral  (from  [10]) 


f°e-cZu2uJ0  (Bu)  du  =  -^e 
Jo  °  2C2 

The  Laplace  transform  of  interest  is 


B2 

c2"2inT.(F!ii)diJ=      X     a  4c2  .  (73> 

2C' 


-a  ■   2(n>  s 
2     a-""    S 


°~"    a<">  (74) 


F7(n)  (s)  = — ^-e 


-'im 


where 

A(n)=2oim2(n)s+ax2  (75) 

[Ref.  4:p.  6].  The  probability  density  function  of  Z-  is 
determined  by  using  the  fact  that  convolution  in  the  time 
domain  translates  to  multiplication  in  the  s  domain: 
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FZi(s)=F^  (S)1F%(S)L-1.  (76) 

The  probability  density  function  of  Z.    is  given  by 

fz(zd)  =Sf"1  [  (Fzu>  (s))1  (FZ(«  (s)  )  L-!]  ,  (77 ) 

where  L  is  the  number  of  hops  per  symbol,  1  is  the  number  of 
jammed  hops,  and  ST1  is  the  inverse  Laplace  transform 
operation.  In  general,  (77)  cannot  be  evaluated  in  closed  form 
and  numerical  subroutines  are  used  to  determine  the  inverse 
Laplace  transform  at  each  point. 

For  the  case  of  all  hops  jammed  or  unjammed,  (77)  reduces 
to 

f Xj- ( zifl) ) -ST1  [( Frtf(  s)  )*] 

/••»  „~<*>>  (78) 

=  /  (F7u»  (s)  )LeSZi*  ds. 

Jo        Zj'm 

The  function  in  the  s  domain  in  (78)  is  obtained  by  raising 
(74)  to  the  L  power 


(n) 


7  J  .,.",„.  2<n)c, 

-tal  (S   J*=( , 


2a  z  s+o  z 


[F-w  (s)]L=( ^-   -)Le 

2a- 

im 

Equation  (79)  is  factored  as 


"■im     LS 

La    2o^lo)s+a/  (79) 
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La 


2(n) 


Zim  -j   2(n) 

^°  im 


La^vi/Uo^)  (80) 


2  S+- 


•  (s+ £ -Le      ** 

2o.  2(/,) 

.un 

The  needed  inverse  Laplace  transform  is  obtained  by  combining 
the  transform  property  that 

r[F(s+a)]  =e~azf  (z)  ,  (81) 

the  inverse  Laplace  transform  relation  (from  [11]) 

c£-i[(A)^]=(Jl)~i   (2/5z)f  (82) 

with  the  substitutions 


a= Vr#  (83) 


and 


La-  2(n)o  2 
2)=—-^ -^.  (84) 

After  simplifying,  the  resulting  closed  form  expression  for 
this  special  case  is 

(-^)^Ox2    ";'  (g,Q.a>)        o2  (85) 

^j"  <Zi   ' 2(n) e  JZ.-1<      2{n)   V^"    )i 

where  IL_1  is  the  modified  Bessel  function  of  order  L-l  and 
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2(n) 

i 

2 


Q(n)=L_^rn_  (g6) 


°x 


[Ref.  4:p.  7]. 

C.   NUMERICAL  INTEGRATION  TECHNIQUES 

A  double  numerical  integration  is  needed  to  evaluate  (48), 
which  can  be  very  time  consuming.  By  taking  advantage  of  the 
properties  of  the  Laplace  transform,  the  double  numerical 
integration  is  reduced  to  a  single  numerical  integration.  The 
amount  of  computing  time  needed  is  greatly  reduced  since 


f  Zjfz  {zi)dzi=^1[-Fz  {s)]  (87) 

Jo    *  s      J 


ro    *  s 

and  (76)  can  be  divided  by  s   as  needed  to  carry  out  the  inner 
numerical  integration. 

The  use  of  the  integration  technique  illustrated  in  (87) 
is  the  primary  reason  a  union  bound  solution  is  calculated 
instead  of  an  exact  solution.  As  seen  in  (44),  the  union  bound 
solution  requires  implementation  of  (87)  only  once.  The  exact 
solution,  as  seen  in  (51),  requires  using  (87)  M-l  times 
since,  for  the  DFT  based  receiver,  the  bins  that  ideally 
contain  no  signal  will  in  general  have  different  probability 
density  functions  due  to  Doppler  shift  and  non-uniform 
windowing.  Although  using  (87)  to  integrate  the  probability 
density  function  significantly  reduces  computing  time  compared 
to  direct  numerical  integration,  the  amount  of  computing 
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resources  needed  to  calculate  the  exact  solution  is  much  more 
than  that  needed  to  determine  the  union  bound  solution.  The 
fact  that  only  a  limited  amount  of  these  resources  was 
available  to  generate  results  justified  the  decision  to 
evaluate  the  union  bound  solution. 

The  remaining  integration  is  accomplished  using  Simpson's 
rule,  which  is  chosen  because  of  its  accuracy  as  compared  with 
the  trapezoidal  rule.  Figure  3.1  shows  how  Simpson's  rule  is 
applied  to  integration.  By  dividing  the  z-axis  into  equally 
spaced  elements  and  evaluating  the  function  at  each  of  these 
points  a  close  approximation  of  the  integral  can  be  obtained 
through  the  relation 

f'"f(z)  dz*-^-  (f(z±)  +4f{z2)  +2f(z3)  +4f(z4)  +2f  (z5) 
+  .  .  .  +4f(zN_1)  +f{zN)  ) 

where  an  odd  number  of  points  must  be  used  [Ref.  12 :p.  350]. 
By  choosing  an  endpoint,  varying  the  spacing  between  points, 
and  repeatedly  evaluating  the  integral  until  no  difference  is 
observed  between  successive  iterations,  an  optimal  Az  is 
obtained.  A  close  approximation  to  the  semi-infinite  integral 
is  obtained  by  using  that  value  of  Az  while  increasing  the 
value  of  the  upper  limit  on  the  integral  until  no  difference 
is  detected  between  successive  iterations. 
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Figure  3.1:  Example  of  dividing  a  function  into  egually 
spaced  segments  suitable  for  use  with  Simpson's  rule. 
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IV.   RESULTS 

A.   PREVIOUS  RESULTS 

Studies  concerning  the  effects  of  non-uniform  windowing  on 
a  DFT  based  receiver  in  the  presence  of  a  Doppler  shifted 
signal  have  been  performed  and  documented  in  [  1  ]  and  [2].  Each 
study  draws  conclusions  concerning  the  effectiveness  of  the 
DFT  based  receiver  under  different  circumstances.  The  results 
of  these  studies  are  included  in  order  to  better  understand 
the  conclusions  of  this  thesis. 

Reference  [1]  studies  the  effects  of  non-uniform  windowing 
on  a  DFT  based  receiver  in  a  nonfading  channel  with  a  Doppler 
shifted  signal,  where  the  probability  of  bit  error  as  a 
function  of  the  signal-to-noise  ratio  is  obtained.  It  is 
demonstrated  that  using  a  non-uniform  window  (specifically,  a 
Hamming  window)  results  in  a  degradation  of  performance  as 
compared  with  the  performance  of  a  uniform  window  for  the 
cases  of  little  of  no  Doppler  shift  on  the  incoming  signal. 
However,  performance  enhancement  is  obtained  as  Doppler  shift 
increases.  For  the  case  of  M=8,  a  Doppler  shift  of  £=.33  is 
sufficient  to  warrant  the  use  of  non-uniform  windowing. 
Furthermore,  the  performance  improvement  increases  as  the 
amount  of  Doppler  shift  in  the  signal  increases.  [Ref.  1:  pp. 
2017-2018] 
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The  effects  of  a  Rician  fading  channel  on  the  DFT  based 
receiver  are  investigated  in  [2].  Another  parameter  for 
analysis  is  the  direct-to-diffuse  ratio  discussed  in  the 
previous  chapter.  In  a  manner  similar  to  that  developed  in 
[1],  data  is  generated  relating  the  probability  of  bit  error 
to  the  signal-to-noise  ratio  for  different  direct-to-diffuse 
ratios.  For  large  direct-to-diffuse  ratios  (i.e.,  conditions 
corresponding  to  those  studied  in  [1]),  the  results  correlate 
well.  However,  as  the  direct-to-diffuse  ratio  is  decreased, 
the  non-uniform  window  gives  better  performance  than  a 
rectangular  window  only  when  the  Doppler  shift  decreases.  The 
presence  of  a  strong  diffuse  signal,  as  illustrated  by  a  small 
direct-to-diffuse  ratio,  limits  the  benefits  of  using  a 
rectangular  window. 

B.   EFFECTS  OF  FREQUENCY-HOPPING  AND  PARTIAL-BAND  INTERFERENCE 

Graphs  of  the  worst  case  probability  of  bit  error  as  a 
function  of  signal-to- jamming  noise  ratio  for  fixed  signal-to- 
thermal  noise  ratios  are  constructed  for  one,  two,  three,  and 
four  hops  per  symbol.  These  results  are  obtained  by 
numerically  evaluating  (40),  (48),  and  (49)  for  values  of  y 
ranging  from  1.0  to  0.001  and  retaining  the  worst  case 
performance  for  each  value  of  Eb/N-.  These  results  are  plotted 
in  Figs.  4.1  through  4.4.  The  signal-to-noise  ratio,  direct  to 
diffuse  ratio,  and  the  modulation  order  are  chosen  to  enable 
direct  comparisons  with  the  results  obtained  in  [4]  which 
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Figure  4.1:  Worst  case  probability  of  bit  error  versus 
signal -to- jamming  noise  ratio  for  L=l,  M=4,  direct-to- 
diffuse  ratio=10,  and  SNR=13.35  dB. 
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Figure  4.2:  Worst  case  probability  of  bit  error  versus 
signal-to- jamming  noise  ratio  for  L=2,  M=4 ,  direct-to- 
diffuse  ratio=10,  and  SNR=13.35  dB. 
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Figure  4.3:  Worst  case  probability  of  bit  error  versus 
signal-to- jamming  noise  ratio  for  L=3,  M=4,  direct-to- 
diffuse  ratio=10,  and  SNR=13.35  dB. 
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Figure  4.4:  Worst  case  probability  of  bit  error  versus 
signal-to- jamming  noise  ratio  for  L=4,  M=4,  direct-to- 
diffuse  ratio=10,  and  SNR=13.35  dB. 
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analyzes  a  conventional  (non-DFT)  noise-normalized  MFSK 
receiver  with  quadratic  detection.  The  results  obtained  here 
for  the  case  of  no  Doppler  shift  (£=0)  and  a  rectangular 
window  should  be  similar  to  those  obtained  in  [4]. 

Differences  in  the  results  obtained  in  [4]  and  those 
obtained  in  this  study  are  primarily  due  to  the  facts  that  the 
equivalent  system  bandwidths  used  in  [4]  and  this  thesis  are 
different  and  the  results  in  [4]  are  exact  while  the  results 
in  this  study  represent  a  union  bound.  The  equivalent  system 
bandwidth  used  in  [4]  is  chosen  to  be  minimal  while  the 
equivalent  system  bandwidth  used  in  this  work,  is  chosen  to 
enable  increased  frequency  spacing  between  the  DFT  output  bins 
as  well  as  to  allow  comparisons  to  be  made  with  the  results  in 
[1]  and  [2].  From  this,  it  is  expected  that  the  data  generated 
for  this  thesis  will  predict  slightly  worse  performance  than 
that  predicted  in  [4].  This  is  observed  in  all  cases. 

Just  as  reported  in  [4],  the  performance  of  the  DFT  based 
communications  system  is  enhanced  as  the  number  of  frequency- 
hops  per  symbol  is  increased.  This  is  especially  true  for  the 
cases  of  little  or  no  Doppler  shift.  For  cases  of  large 
Doppler  shift  (i.e.,  £=0.5),  the  performance  of  the  receiver 
is  slightly  poorer  as  the  hop  rate  is  increased  when  non- 
uniform windowing  is  used.  The  error  introduced  by  self- 
induced  crosstalk,  overcomes  the  benefits  usually  associated 
with  fast  frequency-hopping.  This  is  clearly  seen  in  Fig.  4.1 
and  Fig.  4.4  for  the  case  of  a  rectangular  window  and  £=0.5. 
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For  Eb/N.=41.0  dB,  Pb=4.41xl0~2  for  a  hop  rate  of  one  hop  per 
symbol,  and  Pb=6.65xl0~2  for  a  hop  rate  of  four  hops  per 
symbol.  The  use  of  a  Hamming  window  in  this  situation  improves 
receiver  performance  as  the  number  of  hops  per  symbol 
increases,  but  the  benefits  are  small  compared  to  the 
complexities  needed  to  implement  a  fast  frequency-hopped 
receiver. 

For  the  case  of  slow  frequency-hopping,  the  results  in 
Fig.  4.1  approach  those  reported  in  [2]  as  Eb/N.  becomes  very 
large.  This  is  expected  because  the  decision  statistic 
consists  only  of  one  component  due  to  the  fact  that  the 
carrier  frequency  is  changed  once  per  symbol.  Furthermore, 
receiver  performance  degrades  as  either  the  jamming  power  or 
the  Doppler  shift  increases.  As  reported  in  [1]  and  [2],  the 
performance  of  the  Hamming  window  is  worse  than  that  for  a 
rectangular  window  for  small  Doppler  shifts,  but  outperforms 
the  rectangular  window  for  large  Doppler  shifts.  The  value  of 
£  at  which  the  use  of  a  Hamming  window  gives  better 
performance  than  the  use  of  a  rectangular  window  (the 
crossover  fractional  frequency  offset)  is  a  function  of  the 
amount  of  jamming  power  at  the  receiver.  At  Eb/N.=1.0  dB,  for 
example,  with  a  rectangular  window,  Pb=0.413  with  £=0  and 
Pb=0.687  with  £=0.5;  while,  with  a  Hamming  window,  Pb=0.524 
with  £=0  and  Pb=0.668  with  £=0.5.  However,  at  Eb/N.=41.0  dB, 
with  a  rectangular  window,  Pb=4.59xlO"A  with  £=0  and 
Pb=4.41xl0"2  with   £=0.5;   while,   with  a  Hamming  window, 
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Pb=1.26xl0"3  with  £=0  and  Pb=7.87xl0~3  with  £=0.5.  For  large 
Doppler  shifts,  the  use  of  a  Hamming  window  clearly  improves 
receiver  performance  over  a  broad  range  of  Eb/N. . 

This  same  phenomena  exists  for  the  cases  of  two,  three, 
and  four  hops  per  symbol.  In  these  cases,  however,  the  larger 
hop  rates  have  little  effect  on  the  value  of  the  crossover 
fractional  frequency  offset.  This  is  illustrated  in  Figs.  4.5 
through  4.8  which  plot  the  worst  case  probability  of  bit  error 
as  a  function  of  the  fractional  frequency  offset  for 
Eb/N.=16.0  dB  and  L=l,2,3  and  4.  These  results  are  summarized 
in  Fig.  4.9,  which  plots  the  crossover  fractional  frequency 
offset  as  a  function  of  the  hop  rate  for  various  values  of 
Eb/N-.  As  can  be  seen,  the  values  of  the  crossover  fractional 
frequency  offset  remain  in  a  small  band  for  each  value  of 
Eh/N-,  and  increase  as  the  amount  of  jamming  power  at  the 
receiver  increases. 

From  Fig.  4.9,  it  can  be  seen  that  the  crossover 
fractional  frequency  offset  does  not  follow  a  smooth  curve  as 
the  hop  rate  is  varied.  In  viewing  these  curves,  it  must  be 
remembered  that  the  results  in  Figs.  4.1  through  4.4  represent 
the  worst  case  performances  of  union  bound  solutions.  These 
worst  case  results  are  obtained  by  varying  the  value  of  the 
fraction  of  spectrum  jammed  (y)  and  retaining  only  the  worst 
case  performance.  The  value  of  y  giving  worst  case  performance 
at  one  hop  rate  is  not  generally  the  same  value  at  a  different 
hop  rate.  The  fact  that  only  a  representative  number  of 
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discrete  values  of  y   are  used  to  approximate  the  continuous 
range  of  possible  y's  contributes  to  this  situation. 
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Figure  4.5:  Worst  case  probability  of  bit  error  versus 
fractional  frequency  offset  for  the  case  of  L=l,  M=4, 
SNR=13.35  dB,  and  Eb/N.=16.0  dB . 
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Figure  4.6:  Worst  case  probability  of  bit  error  versus 
fractional  frequency  offset  for  the  case  of  L=2,  M=4, 
SNR=13.35  dB,  and  Eb/N.=16.0  dB. 
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Figure  4.7:  Worst  case  probability  of  bit  error  versus 
fractional  frequency  offset  for  the  case  of  L=3,  M=4, 
SNR=13.35  dB,  and  Eb/N.=16.0  dB . 
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Figure  4.8:  Worst  case  probability  of  bit  error  versus 
fractional  frequency  offset  for  the  case  of  L=4,  M=4, 
SNR=13.35  dB,  and  Eb/N.=16.0  dB . 
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Figure  4.9:  Crossover  fractional  frequency  offset  versus 
frequency-hops  per  symbol  for  the  case  of  M=4,  and 
SNR=13.35  dB. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

As  the  frequency-hop  rate  of  the  DFT  based  communications 
system  with  noise-normalization  increases,  receiver 
performance  improves  in  the  absence  of  large  Doppler  shifts. 
However,  when  large  Doppler  shifts  exist,  the  errors 
introduced  as  a  result  of  self-induced  crosstalk  counteract 
any  benefits  gained  from  the  use  of  fast  frequency-hopping. 
The  advantages  gained  by  using  a  non-uniform  window  in  the 
presence  of  large  Doppler  shifts  improve  the  situation,  but 
not  enough  to  warrant  the  added  hardware  complexities  inherent 
to  the  application  of  fast  frequency-hopped  spread  spectrum. 

The  amount  of  Doppler  shift  necessary  to  warrant  the  use 
of  a  non-uniform  window  is  primarily  a  function  of  the  jamming 
noise  power  spectral  density.  Holding  all  else  constant,  as  N- 
increases,  the  amount  of  Doppler  shift  necessary  to  justify 
the  use  of  a  non-uniform  window  increases.  At  very  small 
values  of  Eh/N-  the  use  of  a  non-uniform  window  results  in 
degraded  performance  for  all  but  the  most  extreme  Doppler 
shifts.  The  performance  degradation  resulting  from  using  a 
non-uniform  window  in  a  situation  with  little  Doppler  shift, 
however,  is  small  compared  to  the  performance  enhancement 
gained  by  using  a  non-uniform  window  when  large  Doppler  shifts 
exist.  In  general,  in  the  absence  of  all  information  regarding 
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the  nature  of  the  received  signal,  the  use  of  a  non-uniform 
window  is  preferred  in  all  situations. 

For  all  but  the  most  severe  conditions  of  Doppler  shift 
and  jamming  power,  the  receiver  studied  in  this  thesis 
provides  reliable  communications  through  the  use  of  fast 
frequency-hopped  spread  spectrum.  Communications  reliability 
does  not  come  without  a  price.  The  complexities  involved  in 
implementing  a  fast  frequency-hopped  spread  spectrum  system 
are  enormous  and  are  not  worth  the  benefits  in  all  situations. 
As  an  example,  in  order  to  keep  the  same  number  of  samples  per 
frequency-hop  for  the  DFT,  the  total  symbol  duration  time 
becomes  a  linear  function  of  the  hop  rate.  In  this  thesis  the 
hop  time  is  determined  to  be  equal  to  one  second.  For  the  case 
of  one  frequency-hop  per  symbol,  this  translates  to  a  symbol 
duration  time  of  one  second,  but  for  the  case  of  L  frequency 
hops  per  symbol,  this  translates  to  a  symbol  duration  time  of 
L  seconds.  Thus,  a  faster  frequency-hop  rate  translates  to 
slower,  but  in  many  situations,  more  reliable  communications. 

Despite  its  drawbacks,  the  use  of  the  DFT  based  receiver 
may  be  superior  in  situations  where  the  content  of  the  message 
is  more  valuable  than  the  speed  of  delivery.  By  employing  the 
DFT  based  receiver  with  fast-frequency  hopping,  a  system  can 
be  developed  which  has  the  potential  to  deliver  reliable 
communications  with  considerable  anti-jam  protection. 
Situations  relying  on  this  type  of  communications  protection 
are  common  in  the  military  environment. 
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In  order  to  gain  deeper  insight  into  the  relationships 
between  the  fractional  frequency  offset,  the  modulation  order, 
the  hop  rate,  and  the  signal  energy- to- jamming  noise  power 
spectral  density  ratio,  further  study  should  be  conducted. 
This  study  should  emphasize  the  changes  in  the  crossover 
fractional  frequency  offset  observed  while  varying  the  key 
parameters  listed  above.  Also,  the  calculation  of  the  exact 
solution  is  possible,  and  should  be  utilized  in  any  further 
study.  Due  to  limited  computational  resources,  this  was  not 
possible  in  this  thesis. 

Finally,  the  work,  presented  here  is  theoretical  in  nature 
and  should  be  tested  with  actual  hardware.  The  assumptions 
made  in  this  analysis  were  selected  to  ease  the  burden  of 
computational  effort  needed  to  solve  the  problem.  They  may  not 
correspond  to  actual  conditions  needed  for  optimal 
communications  in  a  real  world  system.  Only  building  this 
receiver  and  testing  it  under  realistic  conditions  will  prove 
or  disprove  the  usefulness  of  the  DFT  based,  noise-normalized 
receiver. 
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